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ABSTRACT 

The radio counterparts to the 15-/im sources in the European Large Area ISO Survey 
southern fields are identified in 1.4-GHz maps down to ~ 80 /zJy. The radio - MIR 
correlation is investigated and derived for the first time at these flux densities for a 
sample of this size. Our results show that radio and MIR luminosities correlate almost 
as well as radio and FIR, at least up to z ~ 0.6. Using the derived relation and 
its spread together with the observed 15-/zm counts, we have estimated the expected 
contribution of the 15-/xm extragalactic populations to the radio source counts and the 
role of MIR starburst galaxies in the well known 1.4-GHz source excess observed at 
sub-mjy levels. Our analysis demonstrates that IR emitting starburst galaxies do not 
contribute significantly to the 1.4-GHz counts for strong sources, but start to become 
a significant fraction of the radio source population at flux densities < 0.5 — 0.8 mJy. 
They are expected to be responsible for more than 60% of the observed radio counts 
at < 0.05 mJy. These results are in agreement with the existing results on optical 
identifications of faint radio sources. 

Key words: galaxies: evolution - galaxies: starburst - cosmology: observations - 
infrared: galaxies - radio continuum: galaxies. 



1 INTRODUCTION 

With a thousand times better sensitivity than the IRAS 12-fim 
data, the LW3 observations (in the 12 — 18 (im waveband, centred 
at A =15 fim) with the I SOC AM camera (Cesarsky et al. 1996) 
on the Infrared Space Observatory (ISO, Kessler et al. 1996) have 
allowed us for the first time to perform sensitive surveys of dis- 
tant infrared galaxies (up to z ~ 1.5) in the Mid-Infrared (MIR) 
band. The extragalactic source counts derived from these 15-/xm 
surveys, including large area shallow surveys like ELAIS (Oliver 
et al. 2000) covering the flux density range 0.5 < S15 pm < 150 
mJy (Lari et al. 2001), and small area deep integrations reaching 



Sis. 



: 0.05 — 0.1 mJy (Elbaz et al. 1999), show a strong depar- 



ture from no evolution predictions at low flux densities (< 1 — 2 
mJy; Elbaz et al. 1999; Gruppioni et al. 2002). According to both 
optical identification works (Aussel et al. 1999; Elbaz et al. 1999; 
Pozzi et al. 2003) and theoretical models (i.e. Franceschini et al. 
2001), the sources responsible for the sharp upturn observed in 
the number counts at faint flux densities are mainly star-forming 
galaxies at moderately high redshifts (0.4< 2 < 1.4). 



A tight correlation between Far-Infrared (FIR) and radio 
continuum for star-forming galaxies is locally well assessed over 
a large range of luminosities, from normal spirals to the more ex- 
treme Ultraluminous Infrared Galaxies (ULIGS: L > 10 12 Lq), as 
shown by several authors (Condon 1992; Cram et al. 1998; Yun, 
Reddy & Condon 2001). So far there have been many attempts 
to explain the tightness of this correlation, whose origin is still 
somewhat unclear. It is generally assumed that massive stars are 
responsible for both the UV photons heating the dust, which re- 
radiates in the infrared band, and the acceleration of relativistic 
electrons, producing the radio continuum, after their explosion 
as supernovae. Such a tight local correlation does not necessarily 
hold in the distant Universe and, even in this case, it would be 
interesting to investigate possible variations with redshift of the 
slope and/or normalization of the correlation. Moreover, it is not 
obvious that the radio continuum should correlate with the MIR 
emission (produced by a mixture of stochastic heating from Poly- 
cyclic Aromatic Hydrocarbons (PAH; i.e. Puget & Leger 1989) 
and thermal emission at high temperature) as well as with the 
FIR (resulting from thermal emission of large dust grains at lower 
temperature). 
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The aim of this paper is to study for the first time the radio- 
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MIR correlation for galaxies at cosmological redshifts through a 
statistically significant sample. The considered sample contains 
65 sources detected at both 15 /im and 1.4 GHz in the ELAIS 
southern fields, all identified with galaxies up to z ~0.8. By in- 
troducing the radio— MIR correlation into the evolutionary model 
fitting the 15-/im cxtragalactic source counts (see Gruppioni et al. 
2002), we have then estimated the expected contribution to the 
radio source counts from the MIR star-forming galaxies, whose 
role and importance relative to that of AGN in the 1.4-GHz ex- 
cess of sub-mjy/^tjy radio sources is still matter of debate. 

The paper is structured as follows. In section 2 we present 
our data sample. In section 3 we derive the radio-MIR correlation 
and investigate its rcdshift dependence. In section 4 wc estimate 
the contribution of infrared star-forming galaxies to the radio 
source counts and discuss the results and implications. In section 
5 we present our conclusions. 

Throughout this paper we will assume Ho = 75 km s _1 
Mpc -1 , U m = 0.3 and U A = 0.7. 



2 THE DATA SAMPLES 

The ELAIS survey at 15 /im (Oliver et al. 2000), performed in 
raster mode with the ISOCAM instrument, covers a total area of 
~ 12 deg 2 divided into 4 main fields and several smaller areas. 
One of the main fields, SI, and one of the smaller areas, S2, are 
located in the southern hemisphere. SI is centred at a(2000) = 
00 h 34 m 44 _4« ) 5(2000) = -43° 28' 12" and covers an area of 
2° X 2°, while S2 is centred at q(2000) = 05 h 02 m 24.5 s , <S(2000) 
= -30° 36' 00" and covers an area of 21' X 21'. The 15-^tm 
data in these fields have been reduced and analysed using the 
LARI technique, especially developed for the detection of faint 
sources (Lari et al. 2001), obtaining two samples at > 5cr con- 
taining 462 sources with 0.5 < Si5 Mm i 150 mjy in SI (available at 
http://www.bo.astro.it/~elais/catalogues/ELAIS_CAM_15micron_ 
_S1.TAB) and 43 sources with Sis Mm >0.4 mjy, in the deeper 
field S2 respectively. The 15-^tm source counts have been derived 
by Gruppioni et al. (2002) using the ~320 cxtragalactic sources 
detected over the SI area. 

The whole SI and S2 areas have been surveyed in the ra- 
dio with the Australia Telescope Compact Array to Si,4GHz — 
0.2 and 0.13 mjy respectively (Gruppioni et al. 1999; Ciliegi 
et al. 2002, in preparation) and in several optical bands: SI 
in the R band (to R ~ 22.5) with the ESO/Danish 1.5-m 
Telescope and S2 in the U, B, R and / bands with the WFI 
at the ESO 2.2-m Telescope and in the K' band with SOFI 
at the ESO NTT. The radio catalogue in SI (available at 
http://www.bo.astro.it/~elais/catalogues/ELAIS_RADIO_Sl.TAB) 
consists of 652 1.4-GHz sources detected at the 5<r level, while 
the S2 radio catalogue consists of 75 sources. 

Spectroscopic observations of the optical counterparts of the 
ISOCAM sources were carried out at the 2dF/AAT and ESO 
Danish 1.5-m, 3.6-m and NTT telescopes. In SI all the sources 
with R<21.0 have been spcctroscopically identified (~210 extra- 
galactic objects; La Franca ct al. 2002, in preparation). In S2 
we have spectroscopic informations for 29 sources brighter than 
R ~ 21 (about 68% of the sample; Pozzi et al 2003) obtained with 
the ESO 3.6-m telescope. 



3 THE RADIO-MIR CORRELATION 

In order to investigate the radio— MIR correlation within our ISO- 
CAM sample, we have first cross-correlated the 15-/im and the 
1.4-GHz catalogues (complete at the 5<r level) in SI and S2, find- 
ing 28 and 13 coincidences respectively within a distance of 5 
arcsec. Note that the SI 15-^tm catalogue considered for this anal- 
ysis is the same, conservative, catalogue used to derive the source 



Table 1. 15-^m sources detected at 1.4 GHz 



S (mjy) 


(S) 
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Ndet 
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< 0.5 




10 


(0) 


0.0 


0.5 - 0.9 


0.7 


78 


11 (11) 


14.1 


0.9 - 1.6 


1.2 


120 


19 (17) 


15.8 


1.6 - 2.9 


2.2 


87 


28 (27) 


32.2 


2.9 - 5.3 


3.9 


33 


14 (11) 


42.4 


5.3 - 9.5 


7.0 


15 


12 (11) 


80.0 


9.5 - 17.0 


12.7 


7 


6 (5) 


85.7 


17.0 - 30.6 


22.8 


2 


2 (1) 


100.0 


30.6 - 55.1 


41.1 


2 


2 (1) 


100.0 



counts (see Gruppioni et al. 2002), where 35 possibly spurious 
sources (all with S < 1.5 mjy) have been excluded after a visual 
inspection of their pixel history. These sources, detected above 
the 5<r threshold on the maps obtained through a combination 
of several images, are too faint to be distinguished from noise on 
the single pixel histories without uncertainty (see Gruppioni et 
al. 2002 for further explanation). Then, at each ISOCAM posi- 
tion, we have searched for detection in the radio maps down to 
3 cr, finding 53 additional radio identifications within 5 arcsec 
(46 in SI and 7 in S2), for a total of 94 ISO— radio associations 
with Si. 4 GHz __ 3cr. The number of expected spurious coinci- 
dences, on the basis of the density of ISO and radio sources and 
of the adopted maximum distance (5 arcsec) is of the order of 
one. Seventy-three of these sources have spectroscopic data and 
rcdshift measurements in the spectroscopic data sample available 
by October 2002 (additional spectra taken in the October run 
were not used in this work). 

In Table 1 we present the percentage of radio detections as a 
function of the 15-fim flux density. The first two columns give the 
flux density range and the corresponding average flux density re- 
spectively, whilst the following columns give the total number of 
ISOCAM cxtragalactic objects, the number of radio detections, 
the number of non-AGN radio detections(in brackets) and the 
corresponding fraction (i.e. radio detections over total number of 
cxtragalactic sources in that bin). At high flux densities all the 
15-/im cxtragalactic sources have a radio counterpart, while the 
fraction of radio identifications decreases at lower 15-/im flux den- 
sities. This is due to the fact that our radio maps are not deep 
enough to allow the detection of all our 15-^rni sources, especially 
the fainter ones. Figure 1 (left panel) shows the 1.4-GHz lumi- 
nosity versus the 15-fim luminosity for the 65 ISOCAM galax- 
ies with spectroscopic identification detected in the radio band 
(open circles). As radio K-correction we have applied a power 
law with a slope a=0.7, whilst at 15 fim we have applied the K- 
corrections derived by Franceschini et al. (2001) using different 
template Spectral Energy Distributions (SEDs) for the different 
populations (M82 for starbursts and type 2 AGN, M51 for normal 
spirals) modelling the 15-^tm source counts. 

A formal fit to the observed radio-MIR luminosity correla- 
tion for the radio detected galaxies only (dashed line) yields 

log{L 1A ghz/Lq) = (1.09 ± 0.05)log(L ls Mm /L ) (1) 

-(5.91 ±0.54) 

with a dispersion of ~0.27 dex. Slope, normalization and dis- 
persion arc consistent, within errors, with those of the local de- 
termination of the radio— FIR relation found for IRAS galax- 
ies by Yun, Rcddy & Condon (2001; log(L 1A GHz) oc (0.99 ± 
O.Ql)log(L(;o um); u = 0.26). Since our determination does not 
take into account the radio upper limits (about 3/4 of the total), it 
is more representative of the upper envelope (i.e. the brighter ra- 
dio objects) than of the "real" radio-MIR luminosity distribution. 
In order to include the MIR galaxies not detected in the radio, 
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Figure 1. Le/t - Radio (1.4 GHz) versus MIR (15 fim) luminosity for the 15-fim galaxies in the ELAIS southern fields. Objects detected 
in the radio are represented by empty circles (galaxies with spectroscopic redshift) and triangles (without spectroscopic information, to 
which a redshift has been assigned on the basis of an empirical relation between MIR flux density and z: sec text). Arrows are the 3<r 
radio upper limits with measured redshift, while the small dots are the 3<r radio upper limits with no spectroscopic z. The dashed line 
is the best-fitting relation obtained considering only the radio detections with spectroscopic redshift, while the dot-dashed line is the 
best-fitting relation obtained by taking into account also the radio upper limits (and objects without spectroscopic redshift) through 
the ASURV package (Isobc, Fcigclson & Nelson 1986). Right - Radio to MIR luminosity ratio versus MIR luminosity. Symbols, as well 
as the dashed and dot-dashed lines, are the same as in the left panel. The dotted line is the local radio-FIR correlation converted into 
radio— MIR through the empirical relations between Lpjn and L15 found locally by Elbaz et al. (2002). The dot-dot-dot-dashed 
line is the best-fitting relation found by Garrett et al. (2002) for a radio— MIR sample in the HDF-N. 



we have recomputed the radio-MIR luminosity correlation for the 
entire sample of 331 MIR-sclcctcd galaxies by using ASURV (the 
Survival Analysis Package which uses the routines described in 
Isobc, Fcigclson & Nelson 1986 and takes into account also the 
upper or lower limits in a sample). A redshift has been assigned 
to the 192 sources without spectroscopic information by using 
the empirical correlation (and its spread) between the 15-^tm flux 
densities and redshifts found for our 139 spectroscopically identi- 
fied galaxies: 

log(z) = -(0.68±0.04)-(0.37±0.06)«og(5i 5 Mm [mjy])+G(0, a rel )(2) 

where oy e ; (= 0.25) is the la dispersion of the relation and 
G(0,tr re i) is a Gaussian distribution with centre and width 
ay e ;. For the extragalactic sources not detected in the radio we 
have adopted an upper limit to the 1.4-GHz flux density equal to 
their corresponding 3(7 value on radio maps. Under these assump- 
tions, we have re-determined the radio-MIR correlation through 
ASURV, obtaining: 

log(L 1A G hz/Lq) = (1.08 ± 0.04)/o 9 (Li 5 Mm /L ) (3) 

-(6.07 ±0.41) 

with a dispersion of ~0.34 dex. This relation, shown as dot-dashed 
line in the left panel of figure 1, is somewhat more scattered than 
and lies about a factor of 2 below the previous determination. The 
right panel of figure 1 shows the radio to MIR luminosity ratio 
versus MIR luminosity. Our estimate of the "real" correlation 
(dot-dashed line) has a lower normalization than the local radio- 
FIR correlation (dotted line; extrapolated to radio-MIR as de- 
scribed below), which is instead much closer to our determination 
for detections only (dashed line). The local relation corresponds 



to a local value of the "q" parameter equal to 2.34 (defined as q = 
log(L FIR [W]/(3.75 x W 12 [Hz}) x 1/L 1A GHzlWHz- 1 ]), where 
the FIR flux is defined to be 1.26 X 10" 14 (2.585 6 o Mm +Sioo M m) [ w 
m — 2 ]; sec Condon et al. 1991) and is converted to MIR through 
the empirical relations between Lpm and L15 Mm found for local 
galaxies by Elbaz et al. (2002). In the same figure we also show the 
relation derived (using the same K-corrections and cosmological 
parameters considered in our analysis) from the radio and 15-^tm 
data for 19 ISOCAM sources detected in the WSRT deep radio 
survey of the HFD-N region (see Table 1 in Garrett 2002). This 
radio-MIR correlation (shown as dot-dot-dot-dashed line in the 
figure) has a normalization significantly higher (about a factor of 
5) than our best-fitting relation and is also higher than the local 
one. Note, however, that this relation is derived using only the ra- 
dio detections with spectroscopic 2, without taking into account 
the radio upper limits. Since the number of upper limits is similar 
to that of the detections (which constitute ~40% of the 15-^tm 
sample detected at > 5<r level), it is likely that, as we find for our 
data, the "real" correlation would have a somewhat lower nor- 
malization. Despite this, the comparison of our data with those 
in the HDF-N region suggests a possible change in normalization 
of the radio-MIR correlation at the higher redshifts sampled by 
the HDF-N ISO selected galaxies. Alternatively, the two sets of 
data (ELAIS and HDF-N surveys) could be made consistent with 
each other if the ISO or the radio data in the two fields were on 
a different flux scale (i.e. underestimated ISO or overestimated 
radio flux densities in the HDF-N). 

Despite the difference in normalization between our radio- 
infrared correlation and the local one, which seems to suggest a 
change in the radio-infrared correlation with redshift, the impor- 
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Figure 2. Radio/MIR flux density ratios versus redshift for our 
spectroscopically identified galaxies. Diagonal crosses are the me- 
dian values of the ratio in different redshift intervals (z < 0.1, 
0.1 < z < 0.2, 0.2 < z < 0.3 and z > 0.3), plotted in cor- 
respondence of the median redshift. The dashed line shows the 
expected change with z of Si. 4 GHz/Sis due to the differ- 
ence between the average radio and MIR starburst K-corrections 
(normalised to the median value for our data, corresponding to 
Si. 4 GHz/Sl5 fim = 0.15) at z = 0. 



tant result of our analysis is that radio and MIR luminosities for 
galaxies strongly correlate with each other, almost as well as that 
found for the radio and FIR, and at significantly higher redshifts 
than those explored by IRAS. This implies a possible correlation 
between the PAH emission and the radio (and FIR) luminosity. 
The somewhat larger dispersion with respect to that observed for 
the radio-FIR relation is due to the large spread in the mixture of 
PAH and thermal emission at the high temperatures responsible 
for the radiation observed in the MIR band, or to the complicated 
shape of the galaxy SEDs in that wave-band, which introduces un- 
certainties and significant object-to-object variations in the MIR 
K-correction. 

Figure 2 shows the ratio between 1.4-GHz and 15-^tm flux 
densities as a function of redshift for our spectroscopically iden- 
tified sample of galaxies with radio detections. The dashed line 
shows the expected change of this ratio due to the different K- 
corrections (for starburst galaxies) in the two bands. The nor- 
malization of this curve has been chosen so that approximately 
the same number of objects lie above and below the line. This 
corresponds to a value < SuGHz/Smptm >~ 0.15 at z = 0. This 
value is in good agreement with the value obtained by combining 
the well known 1.4-GHz / 60-fim relation - S^o^m — 127 SuGHz 
(Cram et al. 1998) - and the average 15-^tm / 60-^tm ratio found 
by Mazzei et al. (2001) and Xu (2000) - < 5i 5Mm /S 6 o M m >^ 0.05 
(note however above that our radio detections describe the upper 
envelope — i.e. stronger radio sources — of the radio-MIR corre- 
lation, rather than the entire population). Despite the relatively 
large spread in the values of the radio-MIR flux density ratio, no 
obvious trend with redshift is seen in our data for z > 0.07. All 
the objects below this redshift show values of their radio to MIR 
flux density ratios which are significantly lower than the mean. 
One possible reason for this might be that we have missed some 
extended emission in some of these sources. In fact, about half of 
the very low redshift sources have a weak radio emission (detected 
below the 5<r threshold) therefore, because of the low signal-to- 
noise ratio, their flux density may have been underestimated if 
their radio emission is more extended than the beam (~ 15"). In- 
deed, most of these faint radio sources at low redshift are bright 
and extended in both optical and MIR (~ 40" — 1'). 



4 CONTRIBUTION OF MIR GALAXIES TO 
THE RADIO SOURCE COUNTS 

In order to estimate the contribution made to the radio source 
counts from infrared galaxies, we have convolved the correlation 
between radio and 15-^tm luminosities derived in the previous sec- 
tion with the evolutionary model for spiral and starburst galaxies 
which fits the 15-p.m source counts in SI (Gruppioni et al. 2002). 
The intrinsic dispersion adopted for the radio-MIR correlation 
(<r ~ 0.28 dex) has been derived by subtracting from the observed 
value (~0.34 dex) the estimated contributions from the uncertain- 
ties on the radio and 15-/im observed flux densities (~ 0.08 and 
0.10 dex, respectively) and the radio and 15-/im K-corrections 
(~0.02 and 0.15 dex, respectively). The resulting predicted source 
counts for starburst + Seyfert 2, normal spiral and all galaxies, 
are plotted in figure 3 as dashed, dotted and solid lines respec- 
tively. As a consistency check, we have derived the radio counts of 
our ISO selected galaxies directly from the data. Each radio de- 
tected ISO galaxy has been weighted by its radio and MIR spatial 
coverage and the contribution from all the radio detections in ra- 
dio bins have been summed to produce the counts shown as filled 
circles in figure 3. The data points arc in excellent agreement with 
the model predictions; in particular, the star-forming population 
would be responsible for about 40% - 60% of the observed counts 
at Si,4GHz ~ 50 — 100 ^Jy. Therefore, starburst galaxies would 
make up most of the observed radio counts at /tjy level, in agree- 
ment with the results from very deep radio surveys, like that in the 
HDF-N, where ~ 80% of sources detected at S1.4 GHz > 16 /^Jy 
have been identified with starburst galaxies (Richards (2000)). 
Conversely, the fractional contribution of starburst galaxies to the 
radio counts decreases rapidly above SuGHz ~ 0.1 mjy; in good 
agreement with spectroscopic identifications at the sub-mjy level, 
which find ~60-70% of elliptical galaxies and AGN1 among the 
optical counterparts of Si. 4 GHz ~ 0.2 mjy radio sources (Grup- 
pioni, Mignoli & Zamorani 1999). 

Given the differences in the normalization of the relations 
between radio and MIR luminosities derived from different sam- 
ples (see previous section), we have computed the expected radio 
counts by increasing the normalization of our best-fitting correla- 
tion by a factor of 2, in such a way to bring it on the same scale 
as the local correlation. The resulting (total) counts are shown by 
the thin solid line in figure 3. The higher normalization relation 
produces a contribution to the radio counts which is far too high 
with respect to the observed data and exceeds the total radio 
counts at S\.±GHz ~ 0-3 mjy. Of course, this result is produced 
by the combination of the radio-MIR correlation and a model fit- 
ting the source counts. However, the model considered here is the 
only one able to fit the 15-/tm counts in SI, which are lower than 
the other existing ones. Therefore, a radio-MIR correlation with 
a normalization significantly higher than that found in this work 
seems to be inconsistent with the observed faint 1.4-GHz source 
counts. 



5 CONCLUSIONS 

By searching for 3cr detections on deep ATCA 1.4-GHz maps of 
the southern ELAIS fields SI and S2, and cross-correlating these 
with 15-/^m cxtragalactic objects detected respectively by Lari 
et al. (2001) and by Pozzi et al. (2002, in preparation), we have 
obtained a sample of 84 MIR-radio galaxies (65 of which have 
measured redshifts). 

These data have allowed us for the first time to study with 
a statistically significant sample of objects, the radio-MIR corre- 
lation and the radio and infrared properties of galaxies to much 
larger distances and fainter flux densities than previously achieved 
with IRAS. 

The principal results of our analysis are the following: 

© 2002 RAS, MNRAS 000, 1-?? 
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Figure 3. Best-fitting model to our ELAIS 15-^tm extragalactic source counts (in differential form, normalised to a Euclidean non- 
evolving distribution) converted to 1.4 GHz by considering the empirical relation found for our data and its spread. As explained in 
the plot, the thick solid line represents the expected total source counts, while the short-dashed and the dotted lines are respectively 
the modeled contributions of a population of strongly evolving starburst galaxies (plus Scyfcrt 2) and non-evolving spirals. The thin 
solid line is the total contribution from starburst and spiral galaxies obtained by increasing the radio-MIR relation normalization by a 
factor of 2. The filled circles are the counts of the SI 15-/im sources with a > 3cr 1.4-GHz counterpart. The dot-dashed line represents 
the fit of Windhorst, Mathis & Neuschaefer (1990) to the 1.4-GHz counts obtained from 24 different radio surveys. The filled squares 
are the total radio counts in the ELAIS regions (a combination of the SI ATCA data of Gruppioni et al. 1999 with the VLA data in 
the northern ELAIS regions of Ciliegi et al. 1999). The open squares are the 1.4-GHz counts in the HDF-N from Richards (2000). 



1. The radio-MIR correlation for MIR selected galaxies with a 
radio detection in our radio maps is well described by an approxi- 
mately linear relation, with a scatter of ~ 0.27 dcx (similar to that 
found for the local radio-FIR relation), implying that PAH band 
emission correlates with FIR and radio luminosity and that the 
locally determined correlation between radio and infrared emis- 
sion for star-forming galaxies persists to cosmological distances 
(z ~ 0.6). 

2. If wc consider also the radio upper limits we can obtain an 
estimate of the "real" radio-MIR correlation, unbiased by the ra- 
dio non-detections. We still obtain a strong correlation between 
radio and MIR luminosities, but with a factor 2 lower normaliza- 
tion and a larger spread (~ 0.34 dex). The lower normalization 
of our relation corrected for upper limits with respect to the lo- 
cal radio-MIR relation (derived from the radio-FIR through a 
local FIR/MIR average ratio) implies a change in the radio-MIR 
correlation with increasing redshift. 

3. There is no indication of any trend with z in the radio-MIR 
correlation found for our data (apart the K-correction effects), up 
toz~ 0.6. 

4. The contribution of 15-/im galaxies to the radio source counts 
has been computed directly from our data and also by including 
the empirical radio-MIR correlation and its spread into the model 
fitting the MIR extragalactic source counts. Data and model agree 
very well and predict that MIR starburst galaxies should start 
contributing significantly ( > 10%) to the radio counts around 0.5- 
0.8 mjy, with their importance rapidly increasing until they make 
up > 40 — 60% of the observed counts at Si. 4 5 50 — 100 /iJy. 
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